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Current matrix elements and observables for electro- and photo-excitation of baryons 
from the nucleon are studied in a light-front framework. Relativistic effects are exam¬ 
ined by comparison to a nonrelativistic model and can typically be of order 20-25%, 
but can be larger for certain matrix elements, such as radial transitions conventionally 
used to describe the Roper resonance. A systematic study shows that the violation of 
rotational covariance of the baryon transition matrix elements stemming from the use 
of one-body currents is generally small. 

Much of what we know about excited baryon states has grown out of simple non¬ 
relativistic quark models of their structure. These models were originally proposed 
to explain the systematics in the photoconplings of these states, which are extracted 
by partial-wave analysis of single-pion photoproduction experiments. Much more can 
be learned about these states from exclusive electroproduction experiments. Elec¬ 
troproduction experiments measure the Q 2 dependence of these form factors, and so 
simultaneously probe the spatial structure of the excited states and the initial nu¬ 
cleons. Both photoproduction and electroproduction experiments can be extended 
to examine final states other than Ntt, in order to find ‘missing’ states which are 
expected in symmetric quark models of baryons but which do not couple strongly to 
the Ntt channel.0’1 Such experiments are currently being carried out at lower energies 
at MIT/Bates and Mainz. Many experiments to examine these processes up to higher 
energies and Q 2 values will take place at TJNAF. 

It is clear that, once the momentum transfer becomes greater than the mass of 
the constituent quarks, a relativistic treatment of the electromagnetic excitation is 
necessary. However, even at low momentum transfer, the ratio p q /ut q of the average 
quark momentum to its average energy is of order unity, which means that relativis¬ 
tic effects will be significant in any model which describes valence quark degrees of 
freedom. 

The results reported herei make use of light-front Hamiltonian dynamics, 0 in 
which the constituents are treated as particles rather than fields. It shares with 
light-front approaches based upon held theories the property that certain combina¬ 
tions of boosts and rotations are independent of interactions which govern the quark 


dynamics, thus making it possible to perform relatively simple calculations of ma¬ 
trix elements in which composite baryons recoil with large momenta. In addition, we 
make use of a complete orthonormal set of basis states, composed of three constituent 
quarks, which satisfy rotational covariance. Such a basis is the natural starting point 
for dynamical models using the scheme of Bakamjian and Thomas. 0 

A consistent relativistic dynamical treatment of constituent quarks in baryons in¬ 
volves two main parts. First, the three-body relativistic bound-state problem is solved 
for the wavefunctions of baryons with the assumption of three interacting constituent 
quarks. Then these wavefunctions are used to calculate the matrix elements of one-, 
two- and three-body electromagnetic current operators. The conceptual and formal 
background for relativistic, directly interacting quarks is presented in detail in Ref. ill, 
and are outlined briefly below. 

For quantum mechanical systems, relativistic invariance is equivalent to the re¬ 
quirement that there be a consistent set of generators of unitary transformations of 
inhomogeneous Lorentz transformations. For generators of spatial translations (P), 
rotations (J), boosts (K) and time translations (H), that requirement is given by the 
commutation relations. Those relations common to Galilean-invariant systems are 


[J\J k ]=ie jk iJ l - [P M , P v ] = 0 (1) 

[J\ P k } = ie jkl P l ; [J\ K k ] = ie jkl I< 1 (2) 

[K j , P°] = - iP j ; [J j , P°] = 0, (3) 

while those unique to Lorentz-invariant systems are 

[K j , K k ) = -ie jkl J l - [K j , P k ] = iS jk P°. (4) 


An inspection of Eq. (|4|) shows that the Hamiltonian is necessarily linked to at least 
some other generators, e.g., the boosts K. In field theory, the generators are con¬ 
structed via the energy-momentum stress tensor using the exact interacting fields. 
The approach taken here follows that of Bakamjian and Thomas, @ with a direct 
interaction via a mass operator to construct consistent set of generators. 

In light-front dynamics, the generators are reorganized into seven non-interacting 
operators {P + ; Pj_; J 3 ; A" 3 }, and three interacting generators { P ~; Jj_}. The Bakamjian 
construction consists of a mass operator 


M 0 -> M = M 0 + U, 

which determines the interacting generators: 

M 2 + Pj 


P~ = 


P+ ’ 


P+ 


{e 3 x [i(P + - P - )Ej_ - K 3 P J + P ± j 3 + Mj ± } 


(5) 

( 6 ) 


J 


(7) 





( 8 ) 


For a three-quark system, the non-interacting mass operator is 

Mo = Y + 

2=1 

and interactions are added directly to the three-quark system: 

M 0 - M = M 0 (k 1 , k 2 , k 3 ) + U{k u k 2 , k 3 ). (9) 

The Capstick-Isgur interaction,0 consisting of one-gluon exchange plus confining 
terms, satisfies all of the necessary formal requirements for the interaction U. This 
choice violates cluster separability, which is normally a problem for systems of parti¬ 
cles which are individually observable, but not for systems of confined quarks. 

The calculations reported here are described in detail in Ref. i. To calculate the 
current matrix element between initial and final baryon states, we expand in sets of 
free-particle states: 

(M'j;PV|/ + (0)|Mj;P y u) = (2tt) _ 18 y dp' j dp' 2 j c/p 3 J dp 1 J dp 2 j dp 3 Y 

( M'j Py|Pi/iiP 2 /i2P 3 /i 3 ) 
X (Py / lPy / 2Py3l /+ (°)|Pl^lP2^2P3^3) 

x (pi/iip 2 /i 2 P3/i3|Mj; P n). (10) 


The electroweak current operator has a cluster expansion similar to that of its 
nonrelavistic counterpart: 

f'(=) = EW + EW + -- (n) 

3 j<k 

Two-body currents Ijk are required for charge-changing ( e.gi r exchange) and/or non¬ 
local interactions, as they are in the nonrclativistic case, and they are also required 
for covariance of full current. In the front form, one- and two-body currents can 
be grouped separately. We compute only the contributions from one-body matrix 
elements, and assume that the struck quark carries the current of a free Dirac particle: 

(pW + (0)|ph) = 1 (12) 


The baryon state vectors are in turn related to wavefunctions as follows: 


(pi^ip 2j u 2 p 3 /i 3 |Mj; P n) 


d(PnP2,P3) 

0(P,k!,k 2 ) 


(27r) 3 h(p 1 +p 2 + p 3 


X (|hl|h2 1512/^12) (■Sl2hl2|/ i 3| s / i s) 


x (IpfJiphfjLx | Lh l ) (Lh l sij, s I jfJt) 

(k p )y;, M (KA) 4 (*:„/u) 

xdZ 


P) 


( 13 ) 






The quantum numbers of the state vectors correspond to irreducible representations 
of the permutation group. The spins (si 2 , s) can have the values (0,- 7 ), (1, |) and 
( 1 , |), corresponding to quark-spin wavefunctions with mixed symmetry (y p and y A ) 
and total symmetry (x s ), respectively-^ The momenta 


K 

K a 


^(ki-k 2 ); 
-J=(ki + k 2 - 2k 3 ) 


(14) 


preserve the appropriate symmetries under various exchanges of ki, k 2 and k 3 . The 
set of state vectors formed using Eq. ( p~3| ) and Gaussian functions of the momentum 
variables defined in Eq. (|T4|) is complete and orthonormal. Since they are eigenfunc¬ 
tions of the overall spin, they satisfy the relevant rotational covariance properties. 
Any solution to a relativistic model with three constituent quarks can be written as a 
linear combination of these states. Thus, current matrix elements in any such model 
can be expressed in terms of the basis state coefficients and the matrix elements be¬ 
tween basis state vectors. The use of this orthonormal basis allows us to examine the 
transition form factors for many different baryons simultaneously. 

Rotational covariance represents a non-trivial constraint in light-front dynamics. 
It necessitates the existence of two-body current operators because of the interaction 
dependence of the four-vector current. One can test the extent to which rotational 
covariance is violated by constructing a quantity which should vanish under exact 
covariance, and comparing it to non-vanishing physical matrix elements. Helicity 
conservation yields the following constraint: 


Z D ^ABL h )(M'f;P'\V + mMr,P\)Dl(R ch ) = 0, |p'-^| > 2, (15) 


A'A 


where 


Rck = &,( P, BLh = Bc,( P', 


(16) 


IS 


and is a Melosh rotation. For elastic scattering from a nucleon, Eq. 
trivially satisfied. For a transition | —> |, there is a single non-trivial condition, 
while for | 


), there are three. 


The result of combining Eqs. ©-© is a six-dimensional integral over two rel¬ 
ative three momenta. These integrations are performed numerically, as the angular 
integrations cannot be performed analytically. The integration algorithm is a multi¬ 
dimensional quadrature technique recently generalized and extended to higher degree 
and generalized. HI Uncertainties are typically on the order of a few percent for the 
largest matrix elements. The light-quark mass is taken!® to be m u = ntd = 220 
MeV. 

Using the techniques outlined above we can form the light-front current matrix 
elements for nucleon elastic scattering (M N P' ji'\I + {jS)\M N P/x), from Eq.( P~0| ). 


We have evaluated Eq.( |13|) using a simple ground-state harmonic oscillator basis 

state, $ 0 , 0 (k P iK\) = exp j — [k 2 + K'l]/ 2 aj i0 }^ /{^o^ HO ), where the oscillator size 
parameter «ho is taken0-0 to be 0.41 GeV, and using the (Cl) wavefunctions which 
result from the full solution of the relativized model mass operator ,0 expanded up to 
the N — 6 oscillator shell. Eq. ([12]) applies equally well to quark spinor and nucleon 
spinor current matrix elements, so we can extract Fi{Q 2 ) and ^(Q 2 ) for the nucleons 
directly from the above light-front matrix elements. 

Figure |l] compares the proton and neutron Ge and Gm calculated with these two 
wavefunctions, and the modified-dipole fit to the data. Our choice of quark mass for 
the relativistic calculation, while motivated by previous work,00 gives a reasonable fit 
to the nucleon magnetic moments. The single-oscillator relativistic calculation yields 
proton charge and magnetic radii close to those found from the slope near Q 2 =0 of 
the dipole fit to the data. The relativistic calculation using the relativized model 
wavefunctions falls off too slowly with Q 2 , which is due to the larger probability of 
higher-momentum components in these wavefunctions. This confirms the results of 
previous work!! using these wavefunctions, where the nucleon form factors were fit 
by the adoption of relatively soft form factors for the quarks. 

Figure ^ compares the axial-vector form factor Ga(Q 2 ) and Ge(Q 2 1 for the proton 
calculated with the Cl wavefunctions. Relativistic effects are knownQ to reduce the 
axial coupling constant Qa from the static nonrelativistic quark model value of 5/3 
to more like the physical value of 1.25 using simple single-Gaussian wavefunctions; 
using the Cl wavefunctions (]a is reduced flirt, her 12 due to the higher momenta of the 
quarks in these wavefunctions. As expected from the data for Ga(Q 2 ), the axial form 
factor falls with Q 2 more slowly than Ge- 

Figure ^ shows our relativistic results for the A\, A 3 , and C 1 hclicity amplitudes 

22 2 

for electroexcitation of the A| (1232) from nucleon targets using N = 6 Cl wavefunc¬ 
tions for the initial and final momentum-space wavefunctions, compared to relativistic 
results using the single oscillator-basis state above. The parameters ano an d m Ut d 
are the same as above. The relativistic calculation does not solve the problem of 
the long-standing discrepancy between the measured and predicted photocouplings 
(which are essentially transition magnetic moments as the transition is almost purely 
Ml), although the behavior of the single-oscillator relativistic calculation is similar to 
the faster-than-dipole fall off found in the data. Like the nucleon magnetic and axial 
moments, the photocouplings are reduced further by the adoption of the Cl wave- 
functions, and the form factors drop more slowly with Q 2 . A reasonable fit to the Q 2 
dependence of the data is achieved by Cardarclli et alx$ using the Cl wavefunctions 
and soft quark form factors which fit the nucleon form factors. 

We have also plotted the numerical value of the rotational covariance condition 
(multiplied by the normalization factor (Ana / 2 K w for ease of comparison to the 
physical amplitudes), given by the left-hand side of Eq. fll5P , for |/i' — ji\ = 2 . At 




lower values of Q 2 the rotational covariance condition expectation value is a small 

fraction of the transverse helicity amplitudes, but approximately the same size as C 1 

2 

and larger than the value of E2/M1 implied by our Ai and A 3 . 

2 2 

Given the controversy surrounding the nature of the baryon states assigned to 
radial excitations of the nucleon in the nonrelativistic modelS, in Figure [4] we compare 
nonrelativistic and relativistic calculations, for both proton and neutron targets, using 
for the final wavefunction a simple radially excited basis state which can be used to 
represent the P n resonance iV(1440)| + . For the initial state we have used the single 
oscillator-basis ground state wavefunction above. 

There are large relativistic effects, with differences between the relativistic and 
nonrelativistic calculations of factors of three or four. Interestingly, the transverse 
amplitudes also change sign at low Q 2 values approaching the photon point. The 
large amplitudes at moderate Q 2 predicted by the nonrelativistic model (which are 
disfavored by analyses of the available single-pion electroproduction data£3) appear 
to be an artifact of the nonrelativistic approximation. This disagreement, and that of 
the nonrelativistic photocouplings with those extracted from the data for this state,! 
have been taken as evidence that the Roper resonance may not be a simple radial 
excitation of the quark degrees of freedom but may contain excited glue.HE The 
strong sensitivity to relativistic effects demonstrated here suggests that this discrep¬ 
ancy for the Roper resonance amplitudes has a number of possible sources, including 
relativistic effects. 

We also find in the case of proton targets that there is a sizeable C \, reaching a 

2 

1 

value of about 40 — 50 x 10 -3 GeV 2 at Q 2 values between 0.25 and 0.50 GeV 2 , and 
increasing at lower Q 2 values. Correspondingly, there will be a sizeable longitudinal 
excitation amplitude. 

We have also calculated helicity amplitudes for the final state (1535), for 
both proton and neutron targets. Here we use the same single oscillator-basis state 
initial momentum-space wavefunction as above, and final state wavefunctions which 
are made up from momentum-space wavefunctions with one or the other oscillator 
orbit ally-excited. In this case configuration mixing due to the tensor part of the hy- 
perfine interaction is included in the final-state wavefunction. Since the two types of 
orbitally excited states are degenerate in mass before the application of tensor spin- 
spin interactions, they are substantially mixed by them. The results for the helicity 
amplitudes for 7V| (1535) excitation are compared to the corresponding nonrelativis¬ 
tic results in Figure |5|. 

In contrast to the results shown above, here there is reduced sensitivity to rela¬ 
tivistic effects in the results for the transverse amplitudes A 1 / 2 , with the main effect 
being a hardening of the Q 2 behavior of the transition form factor; this is not the 
case for the C\/ 2 amplitudes. For both targets the substantial nonrelativistic C 1/2 


amplitudes are reduced to essentially zero in the relativistic calculation. 

1 Discussion and Summary 

The results outlined above establish that there can be considerable relativistic 
effects at all values of Q 1 2 in the electroexcitation amplitudes of baryon resonances, 
even at Q 2 = 0. In particular, our results show that the Q 2 dependence of the 
nonrelativistic amplitudes is generally modified into one resembling a dipole falloff 
behavior, as has been shown in the case of the nucleon form factors. However, we 
consider it remarkable that relativistic effects account for a large part of discrepancy 
between the nonrelativistic model’s predictions and the physical situation. 

Electroexcitation amplitudes of the Pn Roper resonance iV(1440)| + and 1V(1710)| + 
states, as well as those of the A(1600)§ + , are substantially modified in a relativistic 
calculation. Given the controversial nature of these states,!® we consider this an 
important result. Our results show that relativistic effects tend to reduce the pre¬ 
dicted size of the amplitudes for such states at intermediate and high Q 2 values, in 
keeping with the limited experimental observations for the best known of these states, 
1V(1440)| + . 

We have also found that the rotational covariance violation is a small fraction of 
the larger amplitudes for the Q 2 values considered here. In cases where the dynamics 
causes an amplitude to be intrinsically small, the uncertainty in our results for these 
amplitudes becomes larger. In particular, the calculated ratios E2/M1 and C2/M1 
for the electroexcitation of the A(1232)f + in the absence of configuration mixing 
of D -wave components into the initial and final state wavefunctio nyia are probably 
100% uncertain, and are thus consistent with zero at all Q 2 .0 This may not be the 
case in the presence of such configuration mixing, and we intend to investigate this 
possibility, since A(1232) electroproduction is the subject of current experiments at 
MIT/Bates and several proposed experiments at CEBAF.0 
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by the Department of Energy through Contract No. DE-FC05-85ER250000 (S.C.). 
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Figure 2: Proton axial form factor Ga and Ge calculated with the Cl wavefunctions. 
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Figure 3: Single-basis state and Cl wavefunction relativistic A(1232) electroexcitation 
hclicity amplitudes and rotational covariance condition. 
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Figure 4: Single-basis state relativistic TV(1440) electroexcitation helicity amplitudes, 
and corresponding nonrelativistic Breit-frame helicity amplitudes ( C™/ 2 is zero in the 
nonrelativistic model). 


N(1535), proton target N(1535), neutron target 



Figure 5: Single-basis state relativistic 1V(1535) electroexcitation helicity amplitudes, 
and corresponding nonrelativistic Breit-frame transverse helicity amplitudes. 


















